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We calculate the partition function of the g-state Potts model on arbitrary-length 
cyclic ladder graphs of the square and triangular lattices, with a generalized external 
magnetic field that favors or disfavors a subset of spin values {1,s} with s < q. For 
the case of antiferromagnet spin-spin coupling, these provide exactly solved models 
that exhibit an onset of frustration and competing interactions in the context of a 
novel type of tensor-product S s ®S q - s global symmetry, where S s is the permutation 
group on s objects. 


I. INTRODUCTION 

We consider the g-state Potts model in a generalized external magnetic field that favors 
or disfavors a certain subset of spin values in the interval I s = {l,...,s} with s < q. We 
present exact solutions for the partition function of this model on ladder strips of the square 
and triangular lattices with arbitrary length and periodic (i.e., cyclic) boundary conditions 
in the longitudinal (x) direction and free boundary conditions in the transverse ('//) direction. 
This work extends our previous calculations of the partition function for the cyclic ladder 
ffrapli in the zero-field case in Refs. [1, 0] and for the case of nonzero held and s = 1 in Ref. 

In general, we denote a graph G = (V, E ) by its vertex set V and its edge ( = bond) set E. 
The numbers of vertices, edges, and connected components of G are denoted, respectively, 
by n{G) = n, e(G), and k(G). We note that a triangular-lattice ladder graph of a given 
length may be obtained from the square-lattice ladder graph by adding a diagonal edge in 
each square, connecting, say, the lower left vertex to the upper right vertex in the square. 
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In thermal equilibrium at temperature T, the partition function for the g-state Potts model 
on the graph G in this generalized external magnetic field is given by 


Wi} 


( 1 . 1 ) 


with the Hamiltonian 


'H — J 22 22 22 ^ ai 'P ’ 

(ij) p= 1 i 


( 1 . 2 ) 


where i and j label vertices of G , a* are spin variables with values in the set I q = {1,..., g}, 
(3 = (fc^T) -1 , (ij) denote pairs of adjacent vertices, and J is the spin-spin interaction 
constant. We denote the subset of the q spin values orthogonal to I s as 1^ = {s + 1,..., q}. 
The generalized magnetic field is defined as 


H p = 


H for p E I s 
0 for p G pf 


(1.3) 


Hence, for H > 0 (respectively H < 0), this model favors (resp. disfavors) spin values in 
the interval I s . As is evident in Eqs. (11.11) and (j 1.2 [) . the spin-spin interaction constant J 
and the external magnetic field H always appear in Z multiplied by /3, and, furthermore, in 
exponentiated form, so it is natural to define the following quantities: 


K = (IJ , h = (3H , y — e K , v — y — 1 , w = e h . (1.4) 


The physical ranges of v are v > 0 for the Potts ferromagnet, and — 1 < v < 0 for the Potts 
antiferromagnet. 

There are several motivations for this work. One motivation is the insight that it gives into 
how a new type of global symmetry is manifested in the structure of the partition function. 
If and only if H = 0, then the zero-field Potts model Hamiltonian JL and partition function 
Z are invariant under the global symmetry group in which a,; —> g<Ji V i G V, with g E S q , 
where S q is the symmetric (= permutation) group on q objects. We have used this symmetry 
property to choose, without loss of generality, the set of spins I s to be contiguous. Turning 
on a conventional external magnetic field that favors (or disfavors) one value of the a t has 
the effect of simply reducing the original S q permutation symmetry group to S q -\. (In the 
q — 2 Ising case, this removes any global symmetry, since Si is trivial.) In contrast, with 
the generalized external magnetic held here, setting H ^ 0, reduces the original S q global 
permutation symmetry to the novel type of tensor-product global symmetry group 


£ 


sym 


S S ® Sg- S . 


(1.5) 
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This tensor-product symmetry group in the presence of H ^ 0 has interesting consequences 
for the structure of the partition function, as will be evident from several identities to be 
discussed below. 

A second motivation for the present work is the connection with certain polynomials 
that appear in mathematical graph theory. In the case of zero external magnetic held, the 
partition function of the g-state Potts model is equivalent to a function of central inter¬ 
est in modern mathematical graph theory, namely the Tutte polynomial. If the spin-spin 
exchange constant, J, is negative, i.e., antiferromagnetic, the spin-spin interaction favors 
spin configurations such that spins on adjacent vertices are different. In the limit of this 
Potts antiferromagnet in which the temperature T —> 0 (i.e., v — > —1), the only spin con¬ 
figurations that make a nonzero contribution to the partition function are those such that 
spins on adjacent vertices have different values. It follows that the partition function for 
this zero-temperature g-state Potts antiferromagnet is equal to the chromatic polynomial 
P(G, g), which counts the number of ways of assigning q colors to the vertices of the graph 
G subject to the condition that the colors on adjacent vertices are different. Such a color 
assignment is called a proper g-coloring of (the vertices of) G. Our generalization to a 
nonzero external magnetic field that favors or disfavors a subset of spin values establishes 
a further connection with mathematical graph theory. This connection is embodied in the 
purely graph-theoretic expression that we presented (as Eq. (12.7) in [fj) for the partition 
function. The zero-temperature limit of this partition function in the antiferromagnetic case 
defined what we have called a weighted-set graph coloring problem. Although we will not 
pursue applications here, we have discussed these earlier [4|. For example, the weighted-set 
graph coloring problem could describe the assignment of q frequencies to commercial radio 
broadcasting stations in an area such that adjacent stations must use different frequencies to 
avoid interference and (i) for H > o, the stations prefer frequences in a set I s since reception 
is better for these or (ii) for H < 0, the stations prefer to avoid the frequencies in the set I s 
because of poorer reception. 

A third motivation is that our solutions provide exactly solved examples that exhibit the 
presence or absence of frustration and competing interactions in a controlled, parameter- 
dependent manner. To explain this, we note that for a graph G, the minimum number 
of colors necessary to carry out a proper g-coloring of (the vertices of) G is the chromatic 
number of G, denoted y(G). If the spin-spin exchange interaction, J, is negative, then the 
minimization of the internal energy favors spin configurations in which the spins on adjacent 
vertices have different values. We consider this situation, J < 0, for the remainder of this 
paragraph. Now in the presence of a nonzero external magnetic field H that is positive 
(respectively, negative), the minimization of the internal energy favors spin values that lie 
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ill the interval I s (resp. I^~). The question of whether or not either (or both) of these two 
different situations involves competing interactions and associated frustration depends on 
the value of s or (q — s ) relative to the chromatic number x(G) of G. In the rest of this 
paragraph, we assume that s ^ q since if s = q, then the partition function can be expressed 
completely (see Eq. (12.111) below) in terms of the partition function of a zero-field model 
multiplied by an overall factor. We also assume that s does not have the formal value s = 0, 
since then the H held has no effect (see Eq. (I2.7[) below). Now, if H > 0 and 0 < s < y(G) or 
if H < 0 and 0 < q — s < y(G), then the application of the respective positive or negative H 
involves competing interactions and causes frustration. The reason for this can be explained 
as follows. If 77 > 0 and s < x(G), then the restriction preferred by the magnetic held 
interaction, that the spins should take values in I s , conflicts with the restriction preferred by 
the spin-spin interaction, that adjacent spins should be different, requiring that they take 
on values in the set {1, ...,x(G)}, which is larger than the interval I s . The same statement 
applies for negative H , with the replacement of s by q — s and thus I s by I^~. This frustration 
becomes especially strong as the temperature T —> 0, i.e., K —$■ oo and \h\ —* oo and for 
the infinite strip it leads to a nonanalytic change in the properties of the ground state as a 
function of H/\J\ if J < 0 and s < \{G). A particularly useful aspect of the present work is 
that one can see the origin of this nonanalyticity explicitly from the structure of the exactly 
solved free energy. (For other examples of frustration and competing interactions in spin 
models, see, e.g., [h].) 

More generally, an obvious motivation for exact calculations of partition functions of sta¬ 
tistical mechanical models is the insights that they can give concerning the properties of these 
models that complement those obtained from numerical simulations and series expansions. 

This paper is organized as follows. In Sect. HHwe review basic properties of Z(G, q, s, v, w) 
for an arbitrary graph G. In Sect. IIIII we recall some general structural results for 
Z(G,q, s,v,w) where G is an arbitrarily long cyclic strip graph of the square or triangu¬ 
lar lattice that will be relevant for our specific calculations. In Sects. IIVI and [V] we present 
our calculations of this partition function for arbitrarily long ladder graphs of the square and 
triangular lattices, respectively. In Sect. ED we discuss some physical consequences of these 
calculations. Sect. IVHI contains our conclusions, and several relevant formulas are given in 
the Appendices ITYlfUl 


II. BASIC PROPERTIES 

In this section we will briefly review some basic properties of the partition function that are 
relevant for the present work. Recall that a spanning subgraph G' C G is a graph consisting 
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of the same set of vertices and a subset of the edges of G, i.e., G' = (V, E') with E' c E. It is 
quite useful to have a formula that expresses Z in a purely graph-theoretic manner, instead 
of as a sum over the spin configurations a r . This has been achieved by expressing Z as a sum 
of contributions from the disjoint components Gi comprising spanning subgraphs G' C G. 
The simplest formulation of a g-state Potts model in an external magnetic field H takes this 
field to favor or disfavor one spin value out of the q possible values. In our notation, this 
is the case s — 1 , i.e., the interval I s consists of a single spin value which, without loss of 
generality, can be taken to be a* = 1. For this s = 1 case, a graph-theoretic formula for 
the partition function was given in [gj] (reviewed in [7]). In Refs. j3] and [4, [s| we studied 
this case further and presented a number of properties of the resultant partition function for 
various families of graphs. As discussed in the introduction, a further generalization is to 
take the external magnetic field to favor or disfavor a set of several spin values in an interval 
I s rather than just a single value. For this s > 1 case, we gave a general formula for the 
partition function for an arbitrary graph G (as Eq. (12.7)) in Ref. [4], Refs. ( 9 ], 10] derived 
a number of properties of this partition function and reported exact calculations of it for 
several families of graphs. Let us denote the number of connected components of G as k(G) 
and the connected subgraphs of a spanning subgraph G' as G', i — 1 ,.., fc(G'). Our formula 


is 


Z(G,q,s,v,w) = 


_ KG') 

J2 V< G ' ] Wn(G') 

G'CG i =1 


where 


u m = q — s + sw m = q + s(w m — 1) . 


( 2 , 1 ) 


( 2 . 2 ) 


It follows from Eq. (12 .1 p that Z is a polynomial in the variables q, s, v, and w, and this is 
indicated in the notation Z(G,q, s,v,w). Note that Eqs. (12.ip and (12.2p also allow one to 
generalize q and s from the positive integers to the real numbers. (Indeed, studies of zeros 
of Z(G, q, s, v, w) in q and/or s for fixed v and w require that one generalize q and s further 
to the complex numbers 0, (To| .) Without loss of generality, we restrict ourselves here 
to connected graphs G; however, the spanning subgraphs G' in Eq. (12.lj) may contain (and 
most do contain) more than one connected component. 

The case v = —1 yields a weighted-set chromatic polynomial, denoted Ph(G, q, s,w), 
that describes a graph coloring problem in which one assigns q colors to the vertices of a 
graph such that adjacent vertices have different colors, with a vertex weighting w that either 
disfavors or favors a given subset of s colors contained in the set of q colors (sj-jf]]. In our 
notation, 

Ph(G,q,s,w) = Z(G,q,s,—l,w) . (2.3) 
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As a quantity encoding information about a graph, the weighted-set chromatic polynomial 
Ph(G, q, s, w ) is more powerful than the chromatic polynomial, since it can distinguish be¬ 
tween graphs that are chromatically equivalent, i.e., are different but have the same chromatic 
polynomial [d|- 10|. 

The zero-field g-state Potts model with partition function Z(G,q,v ) is equivalent to a 
function of central importance in mathematical graph theory, the Tutte polynomial [ 111 , 0 ] 
(early reviews include 


13 
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IX 


1 ) 


fc(G')—fc(G) 


(y-i) 


c(G') 


G'CG 


where 


x — 1 + 


(2.4) 


(2.5) 


and y was given in Eq. ('ll.4(1 . so that q = (x— l)(g—1). In Eq. (12.4j) . c(G') denotes the number 
of linearly independent circuits (cycles) in G, satisfying c(G') = e(G') + k(G') — n(G'). Just 
as Ph(G,q,s,w ) generalizes the chromatic polynomial, so also Z(G,q,s,v,w ) generalizes 
the Tutte polynomial. For example, Z(G,q, s, v, w) can distinguish between different graphs 


with the same Tutte polynomnia.1 Urn 


We next review some basic identities and symmetry properties of Z(G,q, s,v,w). From 
the relation w m — 1 = (w — 1) mP with m = n(G[) in conjunction with Eqs. (12.ip and 

(12.21) . it follows that the variable s enters in Z(G,q, s,v,w), and Ph(G,q, s,w) only in the 
combination s(w — 1). Furthermore, Z(G,q,s,v,w ) satisfies the following identities. First, 
if h — 0, i.e., w = 1 , then the partition function reduces to the zero-field expression result, 
i.e., 

Z(G,q,s,v, 1) = Z(G,q,v) . (2.6) 

An equivalent result is obtained if one formally sets s = 0, i.e., the set of spin values that 
interact with the external field is the null set: 


Z(G,q, 0,v,w) = Z(G, q, v) . 


(2.7) 


The equivalence of Eqs. (12.61) and (j2.7[) is in accord with the above-mentioned property that 
Z(G, g, s, v , w) depends on the variables s and w only via the combination s(w — 1), so w — 1 
is equivalent to s = 0. Since the right-hand sides of Eqs. (j2.6[) and (12.71) do not depend on 
s or w, they also imply the partial differential equations 


and 


dPh(G , g, s, 1) 
ds 


( 2 . 8 ) 


dPh{G , g, 0, w ) 
dw 


(2.9) 
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The condition w = 0 means that spin configurations with cq G I s for any i make no contri¬ 
bution to the partition function, leading to the identity 

Z(G,q, s,v,0) — Z(G, q — s,v) . (2.10) 

If s — q, then all of the spin values receive the same weighting due to the generalized external 
magnetic field, so that 

Z(G, q, q, v, w ) = w n Z{G , q, v) . (2.11) 

Finally, from the global tensor-product symmetry (II.5ft of the theory in the presence of the 
generalized magnetic field, there is an identity involving the interchange of the intervals I s 
and I q - S combined with a reversal in the sign of h, namely 

Z(G, q, s, v, w ) = w n Z(G, q,q — s, v, w -1 ) . (2-12) 

Setting v — — 1 in these identities (I2.6p - fl2.12p yields the corresponding relations for 
Ph(G, q, s, w)] for example, Ph(G, q, s, w) = w n Ph(G, q,q — s, ta _1 ), etc. 

III. PARTITION FUNCTION FOR CYCLIC STRIP GRAPHS 

In this section we recall some general structural results from [q| that apply to 
Z(G,q, s,v,w ) for strip graphs of the square or triangular lattice, denoted generically as 
A, of arbitarily great length L x = m vertices and transverse width L y vertices with free 
transverse boundary conditions and periodic (cyclic) longitudinal boundary conditions. We 
denote a strip graph of this type as G = A, L y xm , eye.. The partition function Z(G, q, s, v, w) 
has the general form 

L y nzh{Ly ,d,s) 

Z(k,L y x m,cyc.,q,s,v,w) = ^c (d) (g) ^ ^ w )) m » (3- 1 ) 

d=0 j=l 

where 

8 w («) = c ld \g) = y( 2d ~ 3 )i d ~ i (3.2) 

and q = q — s. The first few ^(q) are c^(g) = 1, 

^(q) = q — l = q — s — 1 , (3.3) 

and 

c^ 2 \q) = q 2 — 3q + 1 = q 2 — 2 qs + s 2 + 3(s — q) + 1 . (3.4) 
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As is evident in the exact solution for the circuit graph in [9] and in our new solutions, 
certain A z,\,L y ,d,j terms occur with multiplicities greater than 1, which are polynomials in s. 
The numbers nzh(L y , d, s) were derived for s = 1 (as Theorem 2.1 with Table 1) in jlij], and 
results for arbitrary s G I s were given in [9|. We have 


^‘ZhkJ-'yi Lyi S) 1 ; (3-5) 

nzhiLy, L y — 1, s) = (s + 1 )L y + ( Ly — l) , (3-6) 

n Z h{L y + 1, 0, s) = (s + 1 )n Z h(L y , 0, s) + n Zh {L y , 1, s) , (3.7) 

and, for 1 < d < L y + 1, 

n Z h(L y + 1 , d, s ) = n Z h(L y , d-l,s) + (s + 2 )n Z h(L y , d, s ) + n Z h(L y , d+l,s) . ( 3 . 8 ) 

For the present case L y = 2, these results yield nzh( 2, 2, s) = 1 as a special case of fj3.5f) . 


n Z h( 2, l,s) = 2s + 3 


and 


? hz/i(2, 0, s) — s 2 + 2s + 2 . 
The total number of A z,sq,L y ,d,j terms is 

Ly 

Nzh,L y ,s = ^ nzh(L y , d, s ) 

d= 0 


' - 5 (?) (t 


If s — 1, then |8j 


For general s [9|, 




s Ly ~ j 


For general s and L y = 2, Eq. (13. lip yields 


N: 


Zh,2,s 


Y] n Z h(,2, d, s) = s 2 + 4s + 6 . 


d =0 


(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 


(3.14) 


A general result is that the unique A z,sq,L y ,d,j for d = L y (for which we thus suppress the 
index j ) is unique, i.e., has multiplicity 1, and is given by 


A Z,sq,Ly,Ly,j — A Z,sq,L y ,Ly — V v 


( 3 . 15 ) 






The s = 1 special case of this was given (as Eq. (4.4)) in [3j. 

The symmetry (12.121) under the replacement of s by q — s and w by w , namely 
Z(G,q,s,v,w) = w n Z(G,q,q — s, U, ru -1 ), has quite interesting consequences for the struc¬ 
ture of this partition function on cyclic lattice strip graphs (of arbitrary length). A basic 
result is that under the replacement s —> q — s, a coefficient c^ d \q) = c^ d \q) changes to 
c^ d \s). Hence, taking into account that the L y x m cyclic strip of the square or triangular 
lattice has n = L x L y = mL y vertices, the symmetry (12.12[) is realized as 

Z( A, L y x m, eye., q, s, v, w ) = w mLy Z{ A, L y x m, eye., q,q — s, v, w -1 ) 

L y nzh{L y ,d,q—s) 

= ^c (d) (s) l wLy ^zA,L y ,dj{q,q - S,n,w _1 )] m . (3.16) 

d =0 j=l 

In order for this identity to hold, there must be an interplay between changes in four different 
quantities: (i) the numbers nzh(L y , d , s ) in Eq. (13.ip which become nzh(L y , d,q — s) in Eq. 
(13.161) : (ii) the related multiplicity factors for various Xz,A,L y ,d,j, which were polynomials in 
s in Eq. (13.11) and become polynomials in q — s in (13.161) : the (iii) the coefficients c^( q ) = 
c^ d \q — s ) in (13.11) , which become c^(s) in (13. 16[) : and (iv) the terms (Xz,A,L y ,d,j) m , which 
also change under the replacements s —» q — s, w —> w^ 1 , and multiplication by w mLy . 

It is natural to analyze how this symmetry (12.121) is realized in the simplest case, namely 
for the circuit graph. Here, nzh( 1,0, s) = s + 1, and nzh{ 1,1, s) = 1. The exact partition 
function for this graph, Z{1D] 1 x m , eye., q, s, v, w) = Z(C m , q, s, v, w) for the circuit graph 
was given in [ 9 ] and is 

s+l 

Z{C m ,q,s,v,w) = Yi X z,Cfi,j(q,s,v,w)] m + c^\q)v m , (3.17) 

3=1 

where the A z,c,d,j terms are as follows. For d — 1, there is a single term 

Xz,c,i,i = v = Xz,c,i ■ (3.18) 


For d — 0 , the Xz,c,o,j consist of the two roots of a quadratic, denoted as Xz.c,o,±, 


where 


A 


1 r 


z,c,o,± 


a±Vr 


A = q + s(w — 1) + v(w + 1) , 
R = A 2 — 4 v(q + v)w , 


(3.19) 
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(3.20) 

(3.21) 




















and, if s > 2, also the terms 


\z,c,o,j — vw for 3 < j < s + 1 . 


(3.22) 


That is, provided that s > 2, the terms A z,c,o,j are all the same for j = 3, ...,s + 1 and 
hence can be considered to be a single A with multiplicity s — 1. Hence, one can re-express 
Z{C m ,q,s,v,w ) in a different but equivalent form (given as Eq. (5.9) in 10]), 


Z{C m ,q,s,v,w) = y^}\z,c,o,j] m + (s - l)(vw) m + (q - s - 1)? 

3 = 1 


(3.23) 


For this circuit graph, the symmetry (j2. 12[) involves the combined effects of changes in 
four quantities: (i) the numbers nzh{ 1,1, s) — 1 and nzh( 1, 0, s) = s + 1, (ii) the multiplicity 
factors for the various terms A z,A,L v ,d,j, hi particular the multiplicity factor s — 1 for a subset 
of the Xz,c,o,j terms; (iii) the coefficient c^(q) = q — s — 1; and (iv) the terms (A z,c,d,j) m 
themselves. Since c^°\q) = 1 and the quadratic roots Xz,c,o,± each occur with multiplicity 
1, it follows that they must transform into w~ x times themselves, i.e., 


s —y q — s 


Xz,c,o,± w l Xz,c,o,± 


(3.24) 


so that with the inclusion of the factor w m in the relevant L y — 1 case of Eq. (13.161) . the 
l(X z ,c,o,+) m + (Az, c, o,-) m ] hi the d = 0 sector is transformed into itself. The other part of 
the d — 0 sector is comprised of the term (s — l)(uia) m coming from the m’th power of 
Az,c,oj = vw with multiplicity s — 1. Under the replacements s —> q — s and w w _1 , this 
term transforms into (q —s— l)(v/w) m , and with inclusion of the w m prefactor, this becomes 
(q — s — l)v m . Under the same replacements s — » q — s and w —> w -1 , the d — 1 term in 
Z(C m ,q,s,v,w), cfX\q)v m = (q — s — l)v m , transforms into (s — l)v m , so with inclusion 
of the w m prefactor, this becomes (s — 1 ){yw) m . Thus, under this replacement, the second 
and third terms in Z(C, q, s, v, w) interchange with each other. This involves what could be 
called a “transmigration” of terms, in which terms from a sector with value of d switch with 
terms in a sector with a different value of d. 

It is also instructive to see how this exact solution satisfies the identity (12.6p for the special 
case of zero field, h — 0, i.e., w — 1. In this case, the s — 1 Xz,c,o,j with 3 < j < s +1 become 
equal to Az,c,i = v. In this w — 1 special case the quadratic roots reduce to q + v and v. 
Thus, (s — 1) —(— 1 = s terms v m associated with the coefficient c^°\q) = 1 become equal to, and 
can be grouped with, the v m term associated with the coefficient c^^(g) = q — s — 1, yielding 
the term v m with coefficient s + (g — s — l) = q — 1 = This shows how the zero-field 
form of the partition function emerges in this limit. In [ 4 J, in the context of the analogous 
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phenomenon for the weighted-set chromatic polynomial Ph(G,q,s,w), we remarked on this 
and denoted it as a “transmigration” process, in which As in a subset of the As associated 
with the coefficient c^ d \q) become equal to a A associated with a coefficient c^ d '\q) with a 
different degree d', so that, taking account of their respective coefficients and multiplicities, 
they can be grouped together. The same phenomenon occurs, in a more complicated manner, 
for Z(sq, 2, q, s, v, w) and Z(tri,2,q,s,v,w). 


IV. PARTITION FUNCTION FOR CYCLIC SQUARE-LATTICE LADDER 

GRAPHS 

Next, we present our exact results for the partition functions of the g-state Potts model 
in a generalized magnetic field on cyclic ladder graphs of the square lattices of length L x = 
m vertices in the longitudinal direction and L y = 2 vertices in the transverse direction, 
Z(sq, 2X171, eye., q, s, v, w ). This graph has n = 2 L x vertices. The partition function has the 
form of (13. ip with nzh{ 2, 2, s) = 1 as a special case of (13. 5p . and with the n^(2, d, s) given in 


Eq. (13.9p for d = 1 and in Eq. (13.1011 for d = 0. As noted above, our new result generalizes 
the zero-field solution given in [1] and the solution for the special case s = 1 given in ^|. 
We calculate the Xz, sq , 2 ,d,jS to be the following: 

• For d = 2, as the L y = 2 special case of (13.15p . there is a single term 


Xz,sq, 2,2,1 = Xz,sq,2,2 ~ V 2 


(4.1) 


For the d = 1 sector, i.e., the \z,sq, 2 ,d,jS with coefficient cd)(g), there are 2s + 3 terms. 
We fold the following results: (i) terms equal to v 2 w, with multiplicity 2(s — 1), which 
thus may be denoted as 

^z,sq,2,i,j = v 2 w for 1 < j < 2(s - 1) , (4.2) 

(ii), next, the terms Xz, S q, 2 ,i,j with j = 2s — 1 and j = 2s, which we denote simply as 
Xz,sq, 2,1,± : 

V[ A±Vr ] , (4.3) 


A 


Z,sq, 2,1,± 


2 L 


satisfying 


Xz,sq, 2,1,± — vXz,C, 1,0, ± J (4-4) 

where Xz,c, i,o,±, A, and R, were given in Eqs. (l3.19K - (l3.2ip . and finally, (iii) the terms 
X sq , 2 ,i,j with 2s + 1 < j < 2s + 3, which are the roots of the cubic equation 


A d -v 


q + s(w — 1) + v(v + 4 + w) 


A 2 
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+ V A 


q(w + 1) + s(w — l)(u + 1) + v{q + 1 + v(w + 1) + 4 w} 


X 


(q + w)(l + v)v 5 w = 0 . 


(4.5) 


The Xz,sq, 2 ,i,± terms in Eq. (14.3ft and the three X z , sq , 2 ,\,j terms that are the roots of 
Eq. (14.51) each have multiplicity 1. 

For d = 0, we have terms Xz }Sq , 2 ,o,j with 1 < j < s 2 + 2s + 2. Since the ordering of the 
index j is a convention, we do not specify it below. These terms and their multiplicities 


are 

( i ) : v(q + v)w with multiplicity 1 , 

(ii) : (vw) 2 with multiplicity s 2 — 3s + 1 , 

two quadratic roots, each with multiplicity s — 1: (iii) 


A 


VVJ 


Z,sq, 2,0,± 


a±Vr 


which satisfy 


Xz,sq,2,0,± — w Xz,sq,2,l,± ~ WV Xz,C,1,0,± , 


(4.6) 

(4.7) 

(4.8) 

(4.9) 


(iv) the three roots of the following cubic equation, each appearing with multiplicity 
s — 1: 


A d — vw 


q + s(w — 1) + v(l + 4 w) + v 2 w 


X 2 


+ (vw) 3 2 q + s(w — l)(u + 1) + v(q + 4 + w) + v 2 (w + 1) 


(■ vw) 5 (q + v)(l + v) = 0 , 


A 


(4.10) 


and (v) the five roots of an algebraic equation of degree 5, each with multiplicity 1. 
This quintic equation is the characteristic polynomial of the transfer matrix Tz )Sq , 2 ,o 
given in Appendix [B] 


In the special case s — 1, Eqs. (j3.9[) and (13. lOjl yield nzh, 2,1 = 5 and nzh, 2,0 = 5, 
respectively. For this special case we have checked that the Xz, S q, 2 ,i,j an d Xz, sq , 2 ,o,j f° r j = 

l, ...,5 agree, respectively, with the eigenvalues that we calculated of the transfer matrices 
Tz,sq, 2 ,i in Eq. (4.18) and T z>sq ^ 0 in Eq. (8.1) of Q. 

It is useful to describe how these Xz, sq , 2 ,d,j terms reduce to the known results for Z(sq, 2 x 

m, q,v) in the zero-field case presented in |1[] (listed, for reference, in Appendix lAl). To 
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distinguish the \ sq , 2 ,d,j terms in the zero-held partition function Z(sq, 2 x m,q,v) from the 
A z, sq , 2 ,d,j terms in the field-dependent partition function Z(sq, 2 x m,q, s,v,w), we suppress 
the subscript Z in the former. Among the Xz, sq , 2 ,d,j with d — 1, first, the \z, sq , 2 ,i,j with 
1 < j < 2(s — 1) obviously reduce to v 2 = X sq , 2 , 2 - The quadratic roots Xz, sq ,i,j with j = 2s — l 
and 2s reduce, respectively as 

Xz, sq , 2 , 1 ,+ = v(q + v ) for w = 1 , (4.11) 

Xz,aq, 2,1 = V 2 for w = 1 . (4.12) 

The roots of the cubic equation (I4.5[) reduce to v 2 and to X sq ^,i,+ in Eq. (IA5I) . each with 
multiplicity 1. 

In the d = 0 sector, the term (i) reduces to v(q + v) with multiplicity 1 and (ii) reduces 
to v 2 with multiplicity s 2 — 3s + 1. The Xz, S q, 2 ,o,± reduce to v(q + v ) and v 2 , each with 
multiplicity s — 1. The roots of the cubic equation (I4.10p reduce to v 2 and A S(?i2 ,i,± hi Eq. 
(IA5j) , each with multiplicity s — 1. Finally, the roots of the quintic equation reduce to v 2 , 
X sq , 2 ,i,± in Eq- (1 A5[l . and A SIJj 2 ,o,± in Eq. (IA6I) . each with multiplicity 1. It is then readily 
conhrmed that in this zero-held case w — 1, our result satisfies the identity 

Z(sq, 2 x m, eye., q , s, v, 1) = Z(sq, 2 x m, eye., q, v ) , (4-13) 

with Z(sq, 2 x m,cyc.,q,v ) calculated in [l|, as given in Eqs. (IA3h - flA5f) . We have also 
conhrmed that our result satishes the various identities given in Sect. HH and that in the 
special case s = 1, it reduces to our previous solution of Z(sq, 2 x m, eye., q, 1, v, w ) given in 


V. PARTITION FUNCTION FOR CYCLIC TRIANGULAR-LATTICE LADDER 

GRAPHS 


In this section we present our exact results for the partition functions of the g-state Potts 
model in a generalized magnetic held on cyclic ladder graphs of the triangular lattice of 
length L x = m vertices in the longitudinal direction, Z(tri, 2 x m,cyc.,q, s,v,w). As was 
the case with the square-lattice ladder strip, the partition function for the triangular-lattice 
ladder strip has the form of (13.ip with %/,(2, 2, s) = 1 as a special case of (13.5p . and with the 


nzh(2,d,s) given in Eq. (13. 9 p for d — 1 and in Eq. (I3.10p for d — 0. Our result generalizes 
the zero-held solution given in Q and the solution for the special case s = 1 given in 
We calculate the A z,tri, 2 ,d,jS to be as follows: 


13 











• For d = 2, there is a single term, 

Xz,tri, 2,2,1 = Xz,tri, 2,2 = , (5-1) 

as the L y = 2 special case of (13.15ft : 

• For the d = 1 sector, i.e., the Xz,tri,2,d,jS with coefficient c^^(g), there are 2s + 3 terms. 
We find the following results: (i) terms equal to v 2 w, with multiplicity 2(s — 1), which 
thus may be denoted as 

Xz,tri,2,l,j = X z ,sq,2,i,j = v 2 w for 1 < j < 2(s - l) , (5.2) 

(ii) next, five terms Xz,tH,2,i,j with 2s — 1 < j < 2s + 3 which are the roots of a quintic 
equation that is the characteristic polynomial of the transfer matrix Tz tri,2,1 given in 
Appendix O 

• For d = 0, there are s 2 + 2s + 2 terms Xz,tri, 2 ,o,j■ Of these, one set coincides with 
the corresponding set for the square-lattice strip, namely (i): (rw) 2 , each occurring 
with multiplicity s 2 — 3s + 1. The remaining terms include (ii) five roots of a quintic 
equation which is the characteristic polynomial of the transfer matrix Tz : tri,2,0a (see 
Appendix [Cj) , each with multiplicity s — 1, and (iii) six terms which are the roots of a 
degree-6 equation that is the characteristic polynomial of the transfer matrix Tz } tri,2,ob 
given in Appendix [UJ each with multiplicity 1. 

The reductions in the case H — 0 (i.e., w = 1) are similar to those for the square-lattice 
strip, so we do not discuss them in detail. 


VI. SOME THERMODYNAMIC PROPERTIES 

A. Thermodynamic Functions 

For a lattice strip graph G m of length L x = m vertices with n = n(G m ) vertices, we 
denote the formal limit of infinite length as 

lim G m = {G} . (6.1) 

ra—>■ 00 

The corresponding dimensionless reduced free energy / in the limit m —> 00 and hence 
n = mLy — > 00 as 

f({G},q,s,v,w) = lim rT l Z(G m ,q,s,v,w) , (6.2) 

n—>00 
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so that the actual Gibbs free energy is Q(T , H) = —ksT f. For the G m = A, 2 xm, eye. cyclic 
strip graphs considered here, in the m —> oo limit, only the Xz.A, 2 .d,j of maximum magnitude 
contributes. This comes from the set of terms with d — 0 [l|, y, |4], [9, 15]. This is in accord 
with the property that / should be independent of the (longitudinal) boundary condition 
and the fact that only the d — 0 sector is present for free longitudinal boundary conditions. 

For the square-lattice ladder strip, among the respective s 2 + 2s + 2 terms with d = 0 
that contribute to Eq. (13.11) . we End that the dominant one is the term of largest magnitude 
among the five roots of the characteristic polynomial of the transfer matrix Tz, sq , 2,0 given 
in Appendix [Bl For the triangular-lattice ladder strip, among the respective s 2 + 2s + 2 
terms with d = 0 that contribute to Eq. (13.11) . we hnd that the dominant one is the term 
of largest magnitude among the six roots of the characteristic polynomial of the transfer 
matrix Tzjh, 2,0 (see Appendix ICj) . 

Thermodynamic functions are obtained from the reduced free energy (j6.2[) in the usual 
way. For first derivatives, the internal energy per vertex is U = — df/d/3 = —Jydf/dy and 
a magnetization per vertex 

Ft f rtf 

(6.3) 


„r df df 

' ~ dh~ W dw 


t is conventional to define a corresponding magnetic order parameter per vertex as (e.g., 

3 ) 

qM — 1 


M = 


Q 


(6.4) 


In the zero-field case, with this definition, in the high-temperature phase, where the global 
Ag-symmetry is realized manifestly, the spins take on random values in the set I q , so M = 1/q 
and hence A4 — 0. In the case of complete ferromagnetic ordering, M = 1, and hence also 
M = 1. 


B. Cases where Frustration and Competing Interactions Occur 

One of the interesting features of our results is that they provide exactly solved examples 
that exhibit the occurrence of frustration and competing interactions as a function of a 
parameter, s, that one can choose. As discussed in the introduction, in the antiferromagnetic 
case with J < 0, if H > 0 (i.e., w > 1) and s < y(G') or if H < 0 (i.e., 0 < w < 1) and 
q — s< x(G), then the imposition of the respective positive or negative H involves competing 
interactions and causes frustration. Given the symmetry (12.12[) . without loss of generality, 
it will suffice to discuss only the case of H > 0, and we shall do so. For our lattice strip 


15 













graphs, 


f 2 if m is even 

x(sq,L y x m,cyc.) = < (6.5) 

13 if m is odd 

f 3 if m — 0 mod 3 

x(tri, L y x m, eye.) = < . (6.6) 

l 4 if m = 1 or m = 2 mod 3 

Since the limit m —> oo can be taken with even m values for the cyclic square-lattice strip 
and with values of m that are equal to 0 mod 3 for the triangular-lattice strip, and since 
only the dominant Xz,A,L y ,o,j is relevant in this limit, it follows that the competition and 
frustration is exhibited in this dominant Xz,A,L y ,o,j for H > 0 and 0 < s < 2, i.e., s = 1 for 
cyclic square-lattice strips and 0 < s < 3, i.e., s = 1 or s = 2, for cyclic triangular-lattice 
strips. In each case, in accordance with our discussion in the introduction, it is assumed 
that in these respective cases, the chromatic number (2 or 3) is no greater than q, so that 
the trivial situation where s = q does not occur. Similarly, there is frustration if H < 0 and 
0 < q — s < 2 for square-lattice strips and 0 < q — s < 3 for triangular-lattice strips. 

Our discussion above is general. An explicit example of this phenomenon is instruc¬ 
tive. The simplest illustration of the frustration is provided by Z for the circuit graph, 
Z(C m ,q,s,v,w). Here, the dominant term is Xz,c,o,+ (given in Eqs. (I3.19p - fl3.2ip ). so 

/({<?}, q, v, s, w) = ln[Az,c,o,+] > ( 6 - 7 ) 


where, as in Eq. (16.11) . C denotes the formal limit of the circuit graph C m as m —> oo. For 
our discussion, it will be useful to have the explicit form of M resulting from Eq. (16.71) . We 
calculate 


M = 


w 

Vr 


2 v(q + v) 


±uyn T u 

s + v - 7 =- , (6.8) 

a + Vr J 

where A and R were given in Eqs. (I3.20p and (13 .21 p . The magnetic order parameter is then 
obtained by substituting this expression for M in Eq. (16.4p . 

In the limit m —> oo, taken on even m, x({C}) = 2. If J < 0 with H > 0, then as T —y 0, 
the spin-spin interaction requires that the spins on adjacent vertices should be different, while 
the magnetic held term biases them to lie in the interval I s . These requirements conflict with 
each other if s — 1, since this is less than %({C'}). In this s = 1 case, as T —> 0, 


Xz,c,o,+ -t - q - 2 + e h |A| + J(q - 2 + e h A K \) 2 + 4 (q - l)e h 


(6.9) 


(where we have used the fact that v — > —1 in this limit). Now, if H/\J\ < 1, i.e., h/\K\ < 1, 
then e h A R \ — y 0 as T — y 0, while e h grows without bound, so Xz,c,o,+ {q — l)e h 

asymptotically, up to exponentially smaller correction terms. If 1 < h/\K\ < 2, then, 
although the e h ~^ term outside the square root and the e 2 ( h ~\ K \) term inside the square root 
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grow, they grow less rapidly than the 4(g — l)e h term inside the square root, so as T —* 0, the 
dominant term in Xz,c,o,+ asymptotically is again y/(q — l)e h . If h/\K\ > 2, then the e 2 (*-—1^1) 
term inside the square root dominates over the 4 (q — l)e h term so, including also the e h ' K ' 
term outside the square root, one finds that, as T —* 0, \z,c,o,+ asymptotically approaches 
e h ~\ K \ with exponentially smaller correction terms. These are the generic intervals in h/\K\. 
One can also trace this demarcation point h/\K\ = 2 separating the two types of behavior in 
the space of complexified parameters, as, e.g. . in another exactly solved spin model exhibiting 
frustration and competing interactions 16]. For completeness, one should also note the 
marginal case where H = 2|J|, i.e., h/\K\ = 2, which forms a set of measure zero in the 
space of the parameters J and H ; in this case, if one expresses \z,c,o,+ i n terms of \K\ alone, 
then as T —> 0, A z,c,o,+ —^ (l/2)e^[l + a/ 4q — 3 ]. 

We now apply this analysis to M, as calculated in Eq. (16.8ft and the resultant A4 from 
Eq. (16.41) . In the case that we are considering, with J < 0 and s — 1, as T —^ 0, if H >2|J|, 
the ferromagnetic tendency due to the external held dominates over the proper (/-coloring 
tendency due to the spin-spin interaction, so 


J < 1, H > 2|J|, s = 1 =>► M = M = 1 . (6.10) 

If, on the other hand, H < 2| J|, then 

J<1, H<2\J\, s = l ^ M = ^, M = 2( g g ~ ^ . (6.11) 

Thus, for J < 0 and s — 1, in the limit T —» 0, as H decreases through the value 2|J|, 
there is a discontinuous decrease in M. from 1 to the value in Eq. (16.111) . For q — 2, this 
latter value is 0; for q — 3, it is 1/4, and as q —> oo, it approaches 1/2 from below. In this 
situation with frustration and competing interactions, there is thus a nonanalytic change in 
the magnetization at T = 0 as a consequence of a change in the ratio H/\J\. Our example 
has the appeal that the physics is manifestly evident because of the explicit exact solution. 


VII. CONCLUSIONS 

In this paper we have presented an exact calculations of the partition functions Z( A, 2 x 
m; q, s, v, w) of the g-state Potts model in a generalized magnetic held, on cyclic ladder strip 
graphs of the square and triangular lattices of arbitrary length. Several interesting analytic 
properties of these results have been elucidated. In particular, we have shown that, for the 
case of antiferromagnet spin-spin coupling, these provide exactly solved models that exhibit 
an onset of frustration and competing interactions, with a resultant nonanalytic change in 
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physical quantities at zero temperature, in the context of a novel type of tensor-product 
S s <g) S q _s global symmetry. 
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Appendix A: Z(sq, 2 x m,cyc.,q,v) 


We review here the result for the partition function Z(sq, 2 x m,cyc.,q,v ) of the cyclic 
square-lattice ladder graph of length L x = m vertices presented in [l[. We include this 
in connection with our discussion in the text showing how our new result for Z(sq, 2 x 
m, eye., q, s, v, w ) reduces to Z(sq , 2xm, eye., q, v ) in the zero-field case w — 1. This partition 
function has the form of 

Ly nz(Ly,d) 

Z(sq, Ly x m, eye., q, v) = ^ c [d) ^ (A sq ,L y ,d,j) m > ( Al ) 

d= o j= 1 

with L y = 2, where n^(2, 0) = 2, nz( 2,1) = 3, nz(2, 2) = 1, and the coefficients are 
given in Eq. (13.21) . so 

c< 0) = 1, = q — 1, c^ 2) = q 2 — 3q + 1 . (A2) 

As is evident in Eq. flAljl . to distinguish the \ sq ,2,d,jS in the zero-field partition function 
Z(sq, 2 x m,q,v ) from the \z,sq, 2 ,d,j in the field-dependent partition function Z(sq, 2 x 
m,q, s,v,w), we suppress the subscript Z in the former. Explicitly, 

Z(sq, 2 x m,cyc.,q,v) = (\ sq ,2,o,+) m + (A sq , 2 ,o-) m 


+ c 


m [(A„, 2 , 1 , 1 )” + (A.,, 2 , 1 ,+)” + (A„,2,1,-)”J + C ( 2 >(A i? , 2 , 2 )” , 

where (in order of decreasing d) X sq ,2,2 = v 2 , 

Xsq,2,l,l = v(q + v) , 


(A3) 


(A4) 


A S g,2,l,(2,3) = A S9 i 2,1,± 
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V 

- 2 

q + v(v + 4) ± (v 4 + 4u 3 + 12u 2 — 2 qv 2 + 4qv + q 2 ) 1 ^ 2 

(A5) 

and 

A S g,2,0 5 (l,2) = ^sg,2,0,=b i \J-^sqdO ) ? 


(A6) 

where 

Asqdo =v 3 + 4v 2 + 3 qv + q 2 


(AT) 

and 

Rsqdo = i’ 6 + 4n 5 — 2q r n 4 — 2g 2 u 3 + 12n 4 + 16gn 3 + 13g 2 n 2 

+ 6g 3 n + q 4 

(AS) 

The reader is referred to [j] for our corresponding solution for the partition function Z(tri , 2x 
m,cyc.,q,v) of the cyclic triangular-lattice ladder graph of arbitrary length. 


Appendix B: T z , sq ,2,o 



Five of the s 2 +2s+2 \ z , sq ,2,o,j terms, each with multiplicity 1, are determined as the roots 
of a quintic equation which is the characteristic polynomial of the transfer matrix Tz tS q,L y ,d 
with L y = 2, d = 0, and the following entries: 


(Tz,sq, 2,o) 1,1 — <T + 3 v(q + v ) 


(Bl) 


{T Z ,sq, 2,o)l,2 = 2 s(q + v)w 


(B2) 


(Tz, sq , 2 ,o) 1,3 — s (v + l)w 2 


(B3) 


(Tz, sq ,2,0)1,4 = s(s - 1 )w 2 


(B4) 


{Tz,s g, 2 , 0 ) 1,5 — {q + 2v)(v + 1) 


(B5) 


{Tz, sg, 2,0)2,1 — Q " + v{2q + v) 


(B6) 


{T z ,sq, 2,0)2,2 = w[ 2 sq + v(q + v)] 


(B7) 


{Tz,sq, 2 , 0 ) 2,3 = (v + l)(n + s)w 2 


(B8) 


(T z ,sq, 2 , 0 ) 2,4 = (v + s)(s - l)w 2 


(B9) 


(2z, 39 , 2 , 0 ) 2 , 5 = (9 + v )( v + 1) 


(BIO) 


(^,39,2,0)3,1 = (2z, 39 , 2 , 0 ) 4,1 = q(q + v) 


(Bll) 


{Tz,sq, 2 , 0 ) 3,2 = (7z,sg,2,o)4,2 = 2g(v + s)w 


(B12) 


(7z, 39 , 2 , 0 ) 3 ,3 = (v + l)(v 2 + 2w + s)w 2 


(B13) 
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(Tz,sq, 2 ,o) 3 , 4 , = (2v + s)(s - l)w 2 (B14) 

{Tz,sq, 2, 0 ) 3 ,5 = (Tz, sq,2, 0 ) 4 ,5 = q(v + 1) (B15) 

{Tz, sq,2,o)i, 3 = + l)(2v + s)w 2 (B16) 

(Tz, 39 , 2 , 0 ) 4,4 = [(s + v) 2 - s - 2y]w 2 (B17) 

(Tz,sq, 2, 0 ) 5,1 = (B18) 

{Tz,sq, 2 , 0 ) 5 ,j = 0 for 2 < j < 4 (B19) 

(Tz, 39 , 2 , 0 ) 5,5 = f 2 ^ + 1) (B20) 


Appendix C: Matrices Tz,tri, 2 ,d for d = 1, d = 0 
1. d = 1 

Five of the s 2 + 2s + 2 Xz,tri, 2 ,i,j terms, each with multiplicity 1, are determined as the 
roots of a quintic equation which is the characteristic polynomial of the matrix Tz,tri,L y ,d 


with L y = 2, d = 1, and the following entries: 

(Tz,tri,2,i)i,i — V Q + 2u 2 (Cl) 

{Tz,tri,2, l)l ,2 = SWV (C2) 

{Tz,tri, 2,l) 1,3 = {Tz,tri,2, 1)2, 3 = ^ (C3) 

(' Tz,tri,2,1 ) 1,4 — {Tz,tri,2,l) 2,4 — 0 (C4) 

(Tz,tri, 2 ,l)l ,5 = {Tz, tri, 2,1 ) 2,5 = v{l + v) (C 5 ) 

{Tz,tri, 2 ,1)2,1 = (C 6 ) 

(Tz,tri, 2 ,1)2 ,2 = WU(U + S) (C7) 

(Tz,tri,2,1)3,1 — (Tz, tri,2,1)4,1 = Vq + 3u 2 + V 3 (C8) 

(Tz, tri, 2 ,1)3 ,2 = (Tz, tri, 2 ,1)4 ,2 = SWV (C9) 

(Tz,tri,2,1)3,3 = vq + 4u 2 + v 3 (CIO) 

(Tz, tri, 2,1)3,4 = SWV (Cll) 

(Tz, tri, 2 ,l) 3,5 ~ {Tz,tri, 2,1 )4,5 — v{l + v){2 + u) (C12) 

(Tz,tri,2,1)4,3 = vq + 3u 2 + v 3 (C13) 
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{Tz,tri, 2 ,l) 4,4 = WV(V + s) (C14) 

(Tz,tri,2,1)5,1 = v 2 q + 3v 3 + v 4 (C15) 

(Tz,tri,2, 1 ) 5,2 = SW *' 2 (C16) 

(Tzjri, 2 , 1 ) 5,3 = 3v 3 + V 4 (C17) 

{Tz,tri, 2 , l) 5,4 = 0 (C18) 

{Tz,tri, 2 , 1 ) 5 ,5 — y2 (l+u)(2 + u) (C19) 


2. d = 0 

Five of the A z,tri,L y ,d,j with L y = 2 and d = 0 are the roots, each with multiplicity s — 1, 
of a quintic equation which is the characteristic polynomial of the transfer matrix Tz t tri, 2 , 0 a- 
This matrix may be obtained from Tz.tr 1,2, 1 by the replacements s —> q — s and w —y w^ 1 
and then multiplication by w 2 . 

Six of the Xz,tri,L y ,d,j with L y = 2 and d = 0 are the roots, each with multiplicity 1, of a 
degree-6 equation which is the characteristic polynomial of the matrix Tz,tri,2,ob with entries 


(Tz,tri, 2 ,ob)i,i ~ T + 4 V Q + 5u 2 + v 3 (C20) 

(Tz,tri,2,0b)l,2 = sw(q + V ) (C21) 

(Tz, tri,2, ob) 1,3 = sw(q + 2v) (C22) 

(Tz, tri,2, ob) 1,4 = -sic 2 (l + v) (C23) 

(Tz,tri,2,0b)l,5 = s ( s ~ T)W 2 (C24) 

(Tz,tri,2,0b)l,6 = (1 + V ){Q + 3v + V 2 ) (C25) 

{Tz,tri,,2,ob) 2,1 — Q 2 + 3 vq + 3v 2 + v 3 (C26) 

(' Tz,tri,2,0b)2,2 = w(v + s)(q + v) (C27) 

(Tz,tri, 2 , Ob) 2 , 3 = Sw(q + v) (C28) 

(Tz,tri,2,Qb)2,4 = VJ 2 (1 + v)(v + s) (C29) 

(Tz,tri,2,Ob) 2,5 = U) 2 (s — l)(u + s) (C30) 

(Tz,tri, 2, Ob) 2,6 = (1 + v)(q + 2v + n 2 ) (C31) 

{Tz,tri, 2 , 06 ) 3,1 — T + + n 2 (C32) 
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{T z ,tri,2,0b)z,2 = wq(v + s) (C33) 

(Tz,tri, 2 ,0fc)3, 3 = w ( v + S )(Q + v) (C34) 

{Tz,tri, 2 , Ob) 3,4 = U> 1 2 (1 + l>)(s + 2v + V 2 ) (C35) 

(Tz,tri,2,Ob) 3,5 = w 2 (2vs -2 v~s + s' 2 ) (C36) 

{Tz,tri,2fib)z,0 = (1 + 1 >)(? + v) (C37) 

(Tz,tri, 2 , 06 ) 4,1 = (Tz,tri, 2, O&Kl = Q(Q + (C38) 

(Tz,tri,2,ob)&,2 = wq(s + 2i> + v 2 ) (C39) 

(Tz,tri, 2 , Ob) 4,3 = (T Z , t.ri, 2 , Ob) 5,3 = Wg(s + v) (C40) 

(Tz,tri, 2 , 0 ^) 4 ,4 = W 2 (l + v)(<S 4“ + 3i> 2 + v 3 4 5 6 ) (C41) 

(Tz,tri, 2 , 0 b) 4,5 = W 2 (s — l)(s + 3v + l’ 2 ) (C42) 

(Tz,tri,2,Ob) 4,6 = (Tz,tri,2,0b)b,6 = Q( 1 + v ) (C43) 

(Tz,tri, 2, ob) 5,2 = + 2-y) (C44) 

(Tz,tri,2,0b)b,4 = W 2 (l + v)(s + 3 ^ + V 2 ) (C45) 

(Tz,tri, 2 , 0 b)b,b = U) 2 (s 2 — S + 3^5 — 3^ + V 2 ) (C46) 

{Tz,tri,2, Ob) 6,1 = v 2 (q + 3v + t? 2 ) (C47) 

{Tz,tri, 2, Ob) 6,2 = Z,tri, 2 , 0 b) 6,4 = (Tz,tri, 2 , Ob) 6,5 = 0 (C48) 

{Tz,tri, 2, Ob) 6,3 = SWV 2 (C49) 

(Tz,t.ri,2, Ob) 6,6 = v 2 (1+v)(2 + v) (C50) 
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